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VIBRATIONS OF A HOLLOW ELASTIC CYLINDER BONDED TO 


A THIN CASING OF A DIFFERENT MATERIAL 
By Vasant S. Kelkar* 

Ames Research Center 


SUMMARY 


Exact solutions are obtained to determine the natural frequencies and 
mode shapes of a thin cylindrical shell supported by a hollow core of a dif- 
ferent material- Materials for both shell and core are assumed to be homo- 
geneous, isotropic, and linearly elastic- A perfect bond is assumed at the 
junction of the shell and the core- The composite cylinder is free from 
stresses at its curved boundaries and is supported by a diaphragm at its flat 
ends. The solutions for the core are based on three-dimensional elasticity 
theory and for the shell on bending theory. Curves are plotted to show the 
variation of the frequency with the variation in circumferential and axial 
wave numbers and in the ratio of inner to outer radii of the core. 


INTRODUCTION 


Problems relating to vibrations of a thick cylinder bonded to a thin 
casing of a different material arise during the flight or transportation and 
handling of solid fuel rocket motors. It is therefore necessary to develop 
solutions which determine the natural frequencies and mode shapes of vibra- 
tions of such composite cylinders. Solid fuel tends to behave like a visco- 
elastic material; however, an analysis based on the assumption of an elastic 
core provides useful results. 

Solutions for some particular cases have been obtained before as cited 
below; however, a completely general solution of the problem using three- 
dimensional elasticity theory for the core and bending theory for the shell 
is difficult to obtain analytically and has not yet been developed. Such a 
general solution, useful in its own right, could also be used to check the 
solutions of approximate methods, such as the finite element method, which 
could then be used to solve problems of more complicated core geometries. Chu 
(ref. 1) gave frequency equations for simple axisymmetric axial shear and 
radial vibrations of composite cylinders, while Achenbach (ref. 2) obtained 
solutions for torsional oscillations . These solutions are relatively simple 
since they involve only one displacement component. Problems involving two of 
the displacement components were also solved; Baltrukonis, Chi, and Gottenberg 
(ref. 3) and Sann and Shaffer (ref. 4 ) obtained frequency equations for plane 
strain vibrations. Neither study considers dis placement variations in the 

*The research was accomplished while the author held a National Academy 
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axial direction which complicate the problem somewhat. Henry and Freudenthal 
(ref. 5) presented solutions for vibrations of a thin shell with a visco- 
elastic core of nonablating boundary. Their solutions are for axisymmetric 
deformations and are based on the use of the correspondence principle. 

For the present analysis, we shall consider a composite circular cylinder 
of length l composed of a thin elastic shell supported by a hollow elastic 
core of a different material (fig. 1). Materials for both shell and core are 
assumed to be homogeneous, isotropic, and linearly elastic, and to be per- 
fectly bonded at their junction. The curved surfaces of the composite cylin- 
der are assumed to be free from stresses. The flat ends of the cylinder are 
assumed to be supported by a diaphragm which prevents displacements in its own 
plane. The solutions presented are completely general and are based on a 
three-dimensional elasticity solution for the core and on bending theory solu- 
tion for the shell. In the analysis which follows, the core and the shell are 
first considered separately. To solve the problem of the composite cylinder, 
equilibrium and compatibility conditions are then satisfied at the junction of 
the shell and the core. The coefficient determinant of the resulting six 
homogeneous equations yields the frequency equation from which numerical val- 
ues of the natural frequency are calculated and plotted against variations in 
axial and circumferential wave numbers for different values of the core thick- 
ness ratio. Frequencies of a composite cylinder with an extremely thin shell 
are compared with those given by Gazis (ref. 6) for a thick cylinder and are 
found to be in good agreement. 


SYMBOLS 


a 

b 


E c> v c 

E s> v s 

e 


ki c 


radius of the middle surface of the shell 
inside radius of the core 


E a 
c 


2(1 + vj 


E s t 

1 - V S * 


elastic constants for the core material 
elastic constants for the shell material 


e x + e cp + € r 


p 2 p 2 - A, 

ar 


2(1 - v c ) 


^P 2 . A 

2v„ a 2 


2 



k s 

l 

m 

n 

P 


^x'^cp^r 

r,cp,x 

t 


-k? and -ki 2 , respectively, when k 2 and k x 2 are less than 
zero 

t 2 

12a 2 

length of cylinder 
circumferential wave number 
axial wave number 
frequency of vibration 

1_ 

pa 

components of applied loading per unit area of shell f s middle 
surface in the x, cp, and r directions 

cylindrical coordinates 

thickness of the shell 

displacements of a point on the middle surface of the shell in 
the x,cp,r directions, respectively 

displacements, respectively, in the axial circumferential and 
radial directions of a point x,cp,r in the core 


W w c displacements of a point on the outer curved surface of the 

core 



A 


njta 

l 


for a cylinder of length 


l 


3 


0- - 2v c )(l + V n ) P 


1 - V, 


C' 1C 

E„ 


c c 

mass density of the material of the core 
mass density of the shell material 


a ,c ,o 
x 9 qr r 


T T T 

rx ; xcp’ rep 


T > T ,0 

rx rep r 


stresses at a point in the core 


stresses at a point on the outer curved surface of the core 


ANALYSIS OF THE CORE 


The core is a hollow circular cylinder for which the governing equations 
for the deformation are (ref. j): 

Equations of motion 


^ T rx , 1 5t *P , a °x , 1 T _ c 
+ + + ^ T r*ir ~ Po. "TT2 


Sr r Sep Sx r rx 
St 


dr r dep dx r r< ^ 


IS. + 1 + + 2 T = 

^ rvn H 


Sa 


Sr^ 


cr - cl 


r + 1 rep + w 'rx + ~r ^ep = 


dt- 

d 2 v 

c 

c dt 2 
d 2 w„ 


dr r dep dx 


dt £ 


( 1 ) 


Kinematic relations 


du. 


€x dx ' 


6 = Xc . 1 ?Ic 

9 r r dtp ^ 

du c dw c 
rx dr dx 


„ = lc + l^c 
*9 dx r dep 


y = 1 + 5lc _ Ic 

/rc P r Sep Sr r 


Sw„ 


r Sr 


( 2 ) 


1 + 



Hooke 1 s law 


where 


x 


cp 


r J 


E, 


1 + v. 


X 


xcp 


rep 


rx j 


E„ 


6c P + 1 - 2 V , 


xcp 


2(1 + V c ) 


rep 


rx 


Su r i Sv P w n Sw~ 
e = £ + £ + e = — — + i + 1 + . > 

X € cp r ax r acp r + a r 


(3) 


Substituting equations (2) into (3) and then introducing the resulting 

expressions for stresses into equations (l), we get three partial differential 

equations in u , v , w , and e. 

c 7 c* 7 c* 7 




1 

_ ^e 

+ V 2 u c = 

2(l + v c ) 


a 2 u c 

( 4 a) 


1 

- 2v 

Sx 

c 

E c 

p c 

at 2 

1 1 Se j 

(v 2 - 

l" 

W + 

2 aw c 

2(1 + v c ) 


^ 2v c 

( 4 b) 

1 - 2v c r acp 

V 

r 2 y 


r 2 Sep 

E c 

p c 

at 2 

<D| 

/ol 

H 

fv 2 - 


> w c " 

2 3v c - 

2(1 + v c ) 


S 2 w 

c 

( 4 c) 

1 - 2v c Sr 

V 

r 2 / 

r 2 Sep 

E c 

p c 

St 2 


where V 2 is the three-dimensional Laplacian operator 

S 2 + 1 _S_ + _1_ b 2 + S 2 
Sr 2 r bv r 2 Sep 2 Sx 2 


The following operations on these equations 


& (■“> + i <*•> + ? & <“ b > + s (1,a) 


give 


5 



(5) 


= Li 2 

at 2 

2 Vq)( 1+ V q) 

1 - v c E c 

We seek a solution of equation (5) which is periodic in cp with a 
period 2jr and periodic in x with a period 2jta/A: 

e = E(r)sin — cos in cpe^^^ (6a) 

a 

where a is a constant which we shall later take as the radius of the middle 
surface of the shell. This gives 


where 


p 2 = 


(i - 


where 


2f§ + 13S+ (y? 

dr 2 r dr \ 

k 2 = p 2 p 2 



The solutions of this equation are 


E(r) — -(l - 2v c )a(ciP m + C 4 R m ) 
where c x and c 4 are arbitrary constants and where 

(1) for k 2 > 0, that is, p 2 > (l/p 2 ) (A 2 /a 2 ) 

I’m - 

and 

«m - 

where J m and Y m are Bessel functions of the first and second kind, 
respectively. 

(2) for k 2 < 0, that is, p 2 < (l/jj 2 ) (A 2 /a 2 ) 

P w = T (Sr) 
m m ' 


(6b) 


(7) 


6 



and 


R* = > k 2 = -k 2 

where I m and are modified Bessel functions of the first and second kind, 
respectively. 

Taking 


u c = f(r)cos — cos mcpe- 5 -^ 
a 

we get from equation (4a) 

+ 1 df + L* . mgN f 1 a E 

dr 2 r dr \ r 2 / 1 - 2v c a 

where 

2 2(l - V c ) 22 A 2 

ki = — up - — 

1 - 2v c a2 

The solutions of the homogeneous part of equation (9) are 


(8) 

(9) 


f(r) = c a% + c 5 S m 

where c 2 and c 5 are arbitrary constants and where 


(10) 


1 - 2 V„ 

(l) for ki 2 > 0, that is, p 2 > 


A fi 


2 „2 


and 


- 2(1 - v c ) H 

= J m( kir ) 

S„ = 


(2) for k x 2 < 0, that is, p 2 < 


1 - 2v 


c 1 A" 


2 n 2 


2(1 - v J V 


Sn kir ) 


and 


S m = V k i r ) > k i 2 = 


-kx c 


7 


When ki 2 > 0* k 2 may be negative or not; hence* to the homogeneous 
solutions of equation (9)> ve can add two different sets of particular solu- 
tions. When ki 2 < 0* k 2 must be less than zero; hence* only one particular 
solution can be added to the homogeneous solutions. All the different 
solutions thus arising can be divided into three cases as follows: 

Case I. kji 2 > k 2 > 0* that is* p 2 /p 0 2 > A 2 > A 2 (l - 2v c )/2(l - v c ) 

Case II. k 2 < 0 < k x 2 * that is* A 2 > p 2 /p Q 2 > A 2 (l - 2v c )/2(l - v c ) 

Case III. 0 > ki 2 > k 2 * that is* p 2 /p 0 2 < A 2 (l - 2v c )/2(l - v c ) < A 2 

where 

P G = l/pa 

Having obtained solutions for e and u c for each of these cases* we can 
solve for the displacements v c and w c from equations (4b) and (4c). This 
procedure has been illustrated for case I in the appendix. Only the solutions 
for the displacements obtained in the three cases are listed here: 

u c = f(r)cos cos mcpe i P t (lla) 

a 

v c = g(r)sin — sin mcpe^^ (lib) 

~ a 

w = h(r)sin 2^ cos mcpe^^ (lie) 

c a 

For a cylinder of length l, A = ma/l (n = 0, 1, . . .,<»). This gives v c , 
w c * and a x as zero at the ends x = 0 and x = l corresponding to a diaphragm 
support. For each f* g* and h* three of the six linearly independent 
solutions are regular and three singular at r = 0. 

Regular solutions 


f(r) 

- AP m } 

Sn 

or 0 



g(r) 

= “ r aP m ’ 

Aim 
a 2 r mi 

or -a 

— Q m 
dr m 

(12) 

h(r) 

= a — P , 
dr m 

- A 1 — Q 
a ki 2 dr 

or a 
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Singular solutions 


f(r) 

~ 4Rjn , 

s m 


or 0 


g(r) 

- - ? aR m > 

A 1 m 
a kx 2 r 

s m 

or " a ^ Sm 

(13) 

h(r) 

= a to ^ , 

A 1 
a kj. 2 

— S m 

dr m 

or a | S m 



where P m , and S m are replaced by the following Bessel functions for 

the various cases • 

C ase Pm ^ m 

I J m (kr) J m (kir) Y m (kr) Y^k^r) 

II I m (kr) Jjkjr) K^fcr) Yjk^r) 

III I m (kr) ImCkir) K^kr) K^r) 

These solutions are the same as those obtained by Gazis (ref. 6) by_using 
scalar and vector potentials. It may be noted that the quantities kr, kr, 
kxr, and k x r are dimensionless, so also are the displacement solutions (ll), 
(12), and (13). We denote solutions (12) by Sj[i(r) [i = 1, 2, 3, j = 1, 2, 

3] and ( 13 ) by &ij(r) [i = 1, 2, 3, j = 4, 5 > 6] • Then the general solution 
for displacements can be written: 

6 

U C 5 ij( r ) c j 

J=1 
6 

V 0 5 a3 (r)Cj 

6 

8 31 (r)C J 
j = l 

where Ci, . . ., c 6 are unknown constants. Corresponding expressions for 
the stresses r rx , T rC p, and a r in the core can be expressed as: 


cos cos m cpe^-^ 


sin sin m cpe ipt 

8L 


(14) 


sin cos m cpe 1 -^ 


9 



where 


6 

T rx = ~2 X a ij( r ) c j C ° S COS m c P elpt 


D 


j=i 

6 


T I^ = ^ tt 2j (r)C j Sin ¥ 31X1 m C f )e±Pt 


D, 


J=1 

6 


a„ = -§■ 7 a 3 .(r)c. sin -^2 CO s m cpe ipt 

J=i 


r _ a 2 L, "3<3^^<3 a 


E a 


D c = 


2(1 + v r ) 


( 15 ) 


The factors aqj(r) can be calculated for the three cases by substituting 
equations (12) and ( 13 ) into equations (2) and (3)* The constants 
ci, • • c 6 can then be determined from stresses and displacements prescribed 
at the inner and outer curved surfaces at the core. This step completes the 
analysis of the core. A few special cases are of interest. 


When m = 0, we get four linearly independent solutions for each f and 
h from (12) and ( 13 ) corresponding to the axisymmetric deformation problem, 
while the two nontrivial solutions obtained for g are for a problem with 
torsional symmetry. 1 

When A = 0, both k 2 and k 2 2 are non-negative. Hence, only case I is 
valid. From solutions (12) and (13), we get four linearly independent solu- 
tions for g and h 1 for the resulting plane strain problem. The two lin- 
early independent solutions for f correspond to pure axial shear vibrations. 


EQUATIONS FOR THE SHELL 


The shell has a midsurface radius a and thickness t which is 
considered to be small compared to a. Using bending theory with Donnell's 
simplifications and neglecting the effect of rotatory inertia, the following 
equations of motion of the shell can be written in terms of its midsurface 
displacements u, v, and w (refs. 8, 9 , and 10). 


1 It is understood that relevant sines and cosines 

solutions (l4) by a simple coordinate transformation. 


are interchanged in 


10 


1 - Vo 

u" + — —£ U* * + 


1 + V 


1 + v 2 

~ v ’’ + v s w + p x §; 


u 1 ' + v’ * + — H v" + w* + p — 

2 2 ^ D g 


v_u' + v’ + w + k (w iv + 2w"‘ * + w**) - p r — 

b b x -n 


a 2 j. d 2 u 


p 

s St 2 


= -p 


a 2 t 6 2 w 
8 D s^ 


(16) 


where 


( )■ . *( ) 

d(x/a) 


and 


( >• = 4 ^ 

dcp 


and p , p , p are loads acting per unit area of the middle surface of the 
shell. X 9 r 


ANALYSIS OF TEE COMPOSITE CYLINDER 


We have a shell of midsurface radius a supported on a core of radius 
a(l - t/2a) which we denote by a Q. We shall satisfy equilibrium conditions 
at the junction of shell and core. Also, for no slip, the displacements of 
the shell and the core must be compatible at their junction. Displacements 
and stresses at the outer surface r = a0 of the core are: 


u c = u c ( x >9>ae) 


Z 6 .(a©)c. cos — cos m cpe^ p ^ 
1 J <3 a 


v c = v c (x,cp,ae) 


w = w (x,cp,ae) 

C ° 


6 

) . ( a0 ) c . sin sin m cpe^"-^ 

/ i J a 


6 


Z S 3 .(a0)c. sin ^ cos m cpe ipt 
0 J a 




( 17 ) 


ii 


I 


6 

7 rx = T rx^^ ,a0 ^ = % Y_, a xj (a0)c O COS C0S m 96 ^ 

j=i 

6 

~rcp = T rC p (x ' cp ' a0) = X a e/ a0)c j Sin IT Sln m Cp ® iPt 

j=i 

6 

o r = o r (x, q>,ae) = % ^ a 3j (ae)c^ sin ^ cos m cpe ipt 

0=1 


(18) 


Displacements u c , v c , w c must be the same as those of a corresponding point 
on the inner surface of the shell. These displacements are related to the 
midsurface displacements u^ v, w of the shell as follows (ref. 10) 


u 


= u + — w 1 
2a 


v 

c 


Gv + 



(19) 


w = w 
c 


We write solutions for u, v, w in the form: 

u = u cos — cos m cpe^ pt 
mn a 

v = t„„ sin M sin a cpe ipt 
mn a 

w = w sin — cos m q)e ipt 
mn a 


( 20 ) 


wherein we may take A = rata /l in order to satisfy diaphragm support condi- 
tions at the ends of a cylinder of finite length l. Substituting 
equations (20) and ( 17 ) into (19) we can express u^, v mn J w mn as 


12 


1 5 ij (a0)c j ~k A 1 V a0)c J 

0=1 0=1 


— / S p . ( a0 ) c . + — — / 5o ^ ( a0 ) c . 
e L 2 J J 2 a e 3 J o 

j=l 0=1 


5 3j( a0 ) C j 


To satisfy equilibrium at the junction of the shell and the core, consider 


the shell to be acted upon on its inner surface by loads 


and -a r 


The moments caused by these forces about axes through the midsurface of the 
shell will be neglected since they will be of rather small magnitude. Hence, 
in equations (l 6) p , p , p r can be directly replaced by -T rx , -T rC p, and 
-a r , respectively. Moreover, since in the derivation of equations (l 6) quan- 
tities of the order t/2a have been neglected compared to 1, we take 0-1 
(ref. 10). When equations (21), (20), and (l8) are substituted into equa- 
tions (l6) , the following three equations result for the constants 

c lj • • v c 6* 


A 2 + 1^ ^-PsX P 2 ) ( 6l j “ k A5s j 


1 + V 


— Am ( i + — m& 3 . 
V 2a 3 J. 


v s A5 3j + 5^ “ij J C j “ ° 

u s - 1 

(22a) 


1 + Vn f 4* 

— Am ( 5, ± - — A5 3i 

2 V 2a 


+ (f + A2 " p s ^ V 2 ) ( 5 aj + k "N 


m5„ . + — a c . =0 
3 J D g j 


(22b) 
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6 


-V S A (63, - — AS ) + m ( s 2 , + — mS 


•J 2a 3 j 


,1 T lUUo J 

-J 2a 3 J 


J=i 


1 + k s (A 4 + 2A 2 m 2 + m 4 ) - p g p £ 

s 


So . + 


^“3,1 

D 3 dJ 


C j = ° 


(22c) 


where 5^ = B^a) and a,^ = a^a). 
inner surface r = b of the core. 


for the constants c^ 


v C 6 


We prescribe zero stresses at the 
This gives us three additional equations 
as follows: 





0 


6 



a £ .( b ) Cj = 0 


I “3 3 (b)cj = 0 

j =1 


(22d) 


(22e) 


(22f) 


Equations (22) together form six homogeneous equations for the six 
constants c^ . . c 6 . They will have a nontrivial solution only if the 

determinant of their coefficients is zero. This gives us the frequency equa- 
tion - a transcendental equation for p, the roots of which can be calculated. 
For each frequency p we can calculate from equations (22) the ratios c 2 /ci, 
C3/ c^ • . Ce/q. The mode shapes are then obtained from equations (±h) , 

( 20 ), and (21). Problems wherein boundary conditions at the inner and outer 
curved surfaces of the composite cylinder are different from those considered 
here can also be solved in an analogous manner. 


NUMERICAL RESULTS 


The results of our analysis are illustrated by numerical values of the 
frequency calculated for case I which covers a wider range of frequency spec- 
trum than case II or III. Numerical values of frequency in these cases can 
be cal ciliated in a similar manner. In case I the dimensionless factors 
ajj(r) have the following form: 


14 



“11 (r) 

“ia(r) 

“is(r) 
“21 (r) 
a 22 (r) 
a 23 (r) 

“ 3 i(r) 

a 32 (r) 

“33(r) 


2 § mAJ m (kr) - 2kaAJ mfl (kr) 


"A 


1 - 


A 


k 1 a a 2 


IT J m (kir) - k 1 aJ TrH . 1 (k 1 r) 


mu' 


mA I J m^ kir ) 


^(m 2 - m) J m (kr) + 2m Mi J mfi (kr) 


s '»<**> - gf W^> 

-2(d^ - m) J§ J m (k!r) + k^a 2 J m ( k x r ) - 2 J^kjr) 


> (23) 


2 ( 1 ^ - m) £- J m (kr) + 2ki- J^^kr) + (A 2 - k^a^jjkr) 


2 {jSL ^ J^xr) 


2A 


ki 2 r 2 


J wl ( k i r ) + 2A J m ( k i r ) 

2 ' - wil r ra ^ 


kir 


2(d^ - m) SJ J (kjr) - 2 j (k x r) 


DH-l' 




The remaining (i = 1, 2, 3> j = 4, 5* 6 ) are obtained by replacing J m 

by Y m and by Y^^ in equations ( 23 ). 

The frequency equation obtained by equating to zero the coefficient 
determinant of equations ( 22 ) was expressed in terms of the dimensionless 
quantity ka. Numerical work was carried out on an IBM 7094 computer. Values 
of v c and v s were fixed at 0.45 and O. 30 , respectively. A relatively weak 
core was considered by taking E c /E s = 10 ” 4 and p c /p s = 0.25. The radius to 
thickness ratio (a/t) of the shell was taken to be 1000. For various values 
of A, m, and a/b, roots ka of the frequency equation were obtained by plot- 
ting the value of the determinant against ka. From the definition of k a 
dimensionless frequency ratio p/p Q is then given by the relation 

(p/p 0 ) 2 = k 2 a 2 + A 2 (24) 

where 

P 0 = l/pa 

In all calculations, the range of values of p/p Q considered was such 
that it satisfied the condition p/p Q > A (or k^a 2 > 0) for case I. The 
results are shown in figures 2 to 6 . 
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For purposes of comparison, numerical values of frequency were also 
obtained for some cases considered by Gazis (ref* 6 ) for the vibrations of a 
thick cylinder. Here, the ratios E c /e s , p c /psj and v c /v s were taken to be 
unity, and a/t was assumed to be 10 4 . With these values it was expected 
that the roots obtained from the frequency equation for our composite cylinder 
would be very close to those given by Gazis for a thick homogeneous cylinder. 
Other parameters assumed to correspond to those used by Gazis were: 

v = v = 0*30 , a/b = 3.0 , A = 1.885 y m = 2 

c s 

The value A = 1.885 corresponds to the cylinder thickness to axial wave 
length ratio of 0.2. Values of p/p , where 


E 


-,1/2 


Pn = 


Acd + v c )J 


were obtained from the frequency equation and compared with those obtained 
from Gazis* curves as follows: 


p/p c 

from frequency 
equation 

— 

O .789 

1.058 

1.415 

1.758 

2.037 

2.257 

p/p c from 
Gasls 1 curves 

0.425 

0.80 

1.07 

1.43 

I .76 

2.06 

2.26 


The first frequency of 0.425 given by Gazis corresponds to a negative 
value of ka and, hence, is out of the range of case I. Since it will be 
obtained from case II or III, it is not given in the comparison above. It 
can be seen that all other roots agree very closely to those given by Gazis 
with a maximum error of 1*375 percent. Similar agreement was also found in 
other cases considered for comparison with Gazis* results. 


DISCUSSION 


For fixed values of the quantities v c , v s and the ratios p c /p s and 
a/t, figures 2 to 6 represent the variation of the dimensionless frequency 
p/p Q for different values of A, m, and the radius ratio a/b. Figures 2 to 
4 show the variation of p/p Q against A for m = 0, 1, and 2, respectively, 
and a/b = 2 . 0 , whereas figures 3 , 5 y and 6 show variation of p/p Q against 
A for a/b = 2, 1.5, and 2.5, respectively. Straight lines p/p 0 = A and 
0.3015 A are drawn on each figure to indicate the domains covered by the 
three cases. However, all the curves are plotted for case I only. Therefore, 
curves are stopped at the dividing line p/p o - I n scmie cases, curves 
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are stopped very close to the dividing line and are extended up to the 
dividing line by dotted lines. 

When m = 0 (fig. 2), the frequency equation degenerates into two 
separate equations, one corresponding to the axisymmetric vibrations (shown 
by solid lines) and the other to pure "torsional" vibrations (shown by dashed 
lines). Furthermore, when A is also equal to zero, the axisymmetric vibra- 
tion frequency equation degenerates further into two equations, one for simple 
radial vibrations and the other for simple axial shear vibrations. When both 
m and A are zero, we expect two of the frequencies to have zero values corre- 
sponding to a rigid body axial translation and rotation of the cylinder. 

Thus, the first two curves in figure 2 go to zero; the dotted curve is for 
torsional motion and the other is expected to be for predominantly 
longitudinal motion. 

Figures 3 and 4 show the variation of frequency with A for m = 1 and 
m = 2. It seems that for higher values of m the frequency variation with 
A is substantially decreased. From figures 3, 5, and 6, one can see the 
effect on the frequency of changing a/b. This change simulates the different 
stages during the burning of the solid fuel. It can be seen that the charac- 
ter of the curves remains the same, however; the magnitude of the frequency 
decreases with increase in a/b ratio. This is as expected since the mass 
increases as a/b increases; however, the stiffness remains substantially the 
same as that of the shell. 


CONCLUDING REMARKS 


Analytical solutions were obtained to determine the natural frequencies 
and mode shapes of a thin cylindrical shell supported by a hollow core of a 
different material. The solutions, obtained by using three-dimensional elas- 
ticity theory for the core and bending theory for the shell, were in closed 
form. Roots of the transcendental frequency equation were obtained numeri- 
cally and Curves were plotted to show a typical variation of the frequency 
with variation in circumferential and axial wave numbers and in the ratio of 
inner to outer radii of the core. The solutions can be used to obtain the 
natural frequencies of cylinders with curved surfaces free from stresses and 
flat ends supported by diaphragms. They can also be used to check frequencies 
obtained by approximate methods, such as finite elements, which, in turn, can 
be used to solve problems of complicated core geometries. 


Ames Research Center 

National Aeronautics and Space Administration 
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APPENDIX A 


SOLUTIONS FOR u,, v Q , AND w Q 


We shall now illustrate the procedure for obtaining the solutions for 
displacements u c , v Q , w c by considering case I in detail. For this case, 
the regular solution for E(r) of equation ( 7 ) becomes 


E(r) = -(1 - 2v)ac 1 J m (kr) 


(Al) 


We shall write only the regular solutions since, for case I, the singular 
solutions can be obtained merely by replacing J by Y. Substituting equa- 
tion (Al) into (9), we get 


+ (*i a - f = CiAJjkr) 
dr 2 r dr \ r 2 / 


(A2) 


The general solution of this equation is given by 


f(r) = c 2 J m (kir) + 


Cl 


k ! 2 - k 2 


AJ m (kr) 


(A3) 


Now from the definition of e 


1 ^Ic 

r Sep 


du r 


W, 


= e - 


dx 


c 

r 


^c 

dr 


(A4) 


Substituting equation (A4) into (4c), we get 


d 2 w 


dw„ 


_c + 3 . 1 

dr 2 r dr ** 


+ -i ^ + 

,2 \ C 


•\2 
O W ( 

dcp" 


d 2 w 


+ 


C + 2 ^c _ 2 e 
.2 r r 


dx^ r dx 

+ — = P c 

1 - 2v„ Sr c 


Se = ft 2(1+ v c ) ^ W c 
E„ St 2 


(A5) 


Substituting for e, u, and w from equations ( 6 a), (Al), (ll), and (A 3 ), we 
get 
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+ 3 to + ( kl 2 

dr 2 r dr \ 



Cl 


2A 2 (1 “ 2v) 

a (_i 2 p 2 


i J m (kr) - 2(1 - 2 v) § J m (kr) 


+ a £ J m (kr) 


+ c 2 ^ J m (kr) 
ar 


(A6) 


Letting H = rh, we can write the left side of equation (A6) as 


+ i as + . s§ j H 

dr 2 r dr V r 2j 


Equation (A6) can then be solved for H; hence. 


h (r) = ? = c x a |L Jjkr) - c 2 ^ ^ A. J m ( kl r) + c 3 f Jj^r) 


h 2 P 2 


(A7) 


Introducing equation (lib) into (A 4 ) one obtains 


g(r) - -=, d - g v)a a Jm(kr) 


d 2 P 2 


+ 


Co A 1 
a ki 2 


2 J (kxr) - c 3 f- J (k^r) 
r x dr in 


(A8) 


Equations (A3), (A 6), and (A8) form three sets of linearly independent solu- 
tions for f, g, and h in case I. Three more sets are obtained by replacing 
the Bessel functions of the first kind by those of the second kind. Simi- 
larly, solutions for the other two cases were also obtained. All these solu- 
tions are summarized (within a multiplicative constant) in equations (12) and 

(13). 
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Figure 6.- Frequency versus A; m = 1, a/b = 2.5. 
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